Effects of Weak Ties on Epidemic Predictability in Community Networks 
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Weak ties play a significant role in the structures and the dynamics of community networks. Based on the 
susceptible-infected model in contact process, we study numerically how weak ties influence the predictability 
of epidemic dynamics. We first investigate the effects of different kinds of weak ties on the variabilities of 
both the arrival time and the prevalence of disease, and find that the bridgeness with small degree can enhance 
the predictability of epidemic spreading. Once weak ties are settled, compared with the variability of arrival 
time, the variability of prevalence displays a diametrically opposed changing trend with both the distance of the 
initial seed to the bridgeness and the degree of the initial seed. More specifically, the further distance and the 
larger degree of the initial seed can induce the better predictability of arrival time and the worse predictability 
of prevalence. Moreover, we discuss the effects of weak tie number on the epidemic variability. As community 
strength becomes very strong, which is caused by the decrease of weak tie number, the epidemic variability 
will change dramatically. Compared with the case of hub seed and random seed, the bridgenss seed can result 
in the worst predictability of arrival time and the best predictability of prevalence. These results show that 
the variability of arrival time always marks a complete reversal trend of that of prevalence, which implies it is 
impossible to predict epidemic spreading in the early stage of outbreaks accurately. 

PACS numbers: 05.40.Fb,05.60Cd,89.75.Hc 



In community networks, the links that connect pairs of 
nodes belonging to different communities are defined as 
weak ties. The weak ties hypothesis, which is first pro- 
posed by Granovetter, is a central concept in social net- 
work analysis. Weak ties not only play a role in effecting 
social cohesion, but also are helpful for stabilizing all com- 
plex systems. Most recent research results showed that 
weak ties have significant impacts on spreading dynamics. 
But until now, no one has given us any study on the effects 
of them on the predictability of epidemic dynamics. In this 
study, we investigate how different kinds of weak ties and 
weak tie number influence the predictability of epidemic 
dynamics in a local community. We show numerically that 
both the degree of bridgeness and the network modularity 
play a significant role in the predictability of the epidemic 
spreading in the local community. More importantly, we 
find that the variability of arrival time always marks a 
complete reversal trend of that of prevalence, which im- 
plies it is impossible to predict epidemic spreading in the 
early stage of outbreaks accurately. This work provides us 
further understanding and new perspective in the effect of 
weak ties on epidemic spreading. 



I. INTRODUCTION 

Community structures at mesoscale level are ubiquitous in 
a variety of real complex systems [1], such as Facebook [2], 
YouTube yfl, and Xiaonei [4]. In general, there are more con- 
nections between members in the same community than be- 
tween members from different communities, where the links 
that connect pairs of nodes belonging to different communi- 



ties are defined as weak ties JD-0]. The weak ties hypothesis, 
which is first proposed by Granovetter JH], is a central concept 
in social network analysis J^, Hol|. W eak ties not only play a 
role in effecting social cohesion lip , but also are helpful for 
stabilizing all complex systems M12I1 . 

Epidemic spreading lfl3l4T6ll . a fundamental dynamical pro- 
cess, is one of the most important subjects in complex network 
theory Il7l - l2lll . Inspired by the significant effects of weak 
ties on network dynamics [ 12], many recent works have con- 
tributed to understanding the interplay between weak ties and 
spreading dynamics in community networks I122l - l26ll . Onnela 
et al. found that weak links can significantly slow diffusion 
process, leading to dynamic trapping of information in com- 
munities 1271 12811 . As weak ties are removed gradually, the 
coverage of information will drop sharply i29ll . In adaptive 
networks, strong communities with weak ties may prevent dis- 
ease propagation ll30l[3lll . 



'Electronic address: tangminghuang521@hotmail.com 



Up to now, almost all studies concentrate only on how weak 
links influence epidemic dynamics in community networks, 
but the predictability of the dynamics is neglected. In order 
to assess the accuracy and the forecasting capabilities of nu- 
merical models, the predictability of outbreaks has been in- 
vestigated in many studies [32- 360 . In view of this point, 
Colizza et al. studied the effect of airline transportation net- 
work on the predictability of epidemic pattern by means of 
the normalized entropy function and found that the het- 
erogeneous weight distribution contributes to enhancing the 
predictability. Crepey et al. found that initial conditions such 
as the degree heterogeneity of the initial seed can induce a 
large variability on the prediction of prevalence Ir38ll . Loecher 
et al. argued that RWC serves as a better index than degree 
to predict the prevalence [39}]. Comparing the scale-free net- 
work (SFN) with community structure with the random SFN, 
the predictability of its global prevalence was found to be 
better [40]. Considering the relative independence of a lo- 
cal community, we investigated the prevalence and its vari- 
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ability in the local community, and found that the extraordi- 
narily large variability in the early stage of outbreaks made 
the prediction of epidemic spreading hard 14111 . In addition, 
we also studied how heterogeneous time delay (HETD) as- 
sociated with geographical distance influences the spreading 
speed and the variability of prevalence. Owing to correla- 
tions between time delay and network hierarchy in HETD, 
epidemic spreading is slowed down obviously and the pre- 
dictability of prevalence is reduced remarkably ll42ll . 

In community networks, as mentioned above, weak ties 
play a very significant role in epidemic dynamics. But un- 
til now, there is no study on the effects of them on the pre- 
dictability of epidemic dynamics. In this paper, we investigate 
how weak ties influence the predictability of epidemic dynam- 
ics in community networks. We show numerically that both 
different kinds and number of weak ties can remarkably in- 
fluence the predictability of the epidemic spreading in a local 
community. More importantly, we find that the variability of 
arrival time always marks a complete reversal trend of that of 
prevalence, which implies it is impossible to predict epidemic 
spreading in the early stage of outbreaks accurately. 

This paper is organized as follows. In Sec. II, we briefly de- 
scribe the dynamical process in community network and intro- 
duce the quantitative measurements of predictability. In Sec. 
Ill, we investigate the effects of different kinds of weak ties 
on the predictability of the dynamics. In Sec. IV, the effects 
of weak tie number on the predictability are analyzed. Finally, 
we draw conclusions in Sec. V. 



H. MODEL INTRODUCTION 
A. Community Network with Degree Heterogeneity 

To investigate the effects of weak ties on the the predictabil- 
ity of epidemic dynamics, we must first identify which links in 
community networks are weak ties. Unfortunately, there is not 
an generally accepted and authenticate community detection 
algorithm 11471 - 14911 , thus it is difficult to identify weak ties ac- 
curately in real networks. We here consider a community net- 
work model comprised of two confined communities A and 
B. Except for community structure, degree heterogeneity is 
another important feature of real community networks IH31 - 
l46ll . In view of this point, we focus on the community net- 
work with degree heterogeneity. To be specific, two indepen- 
dent BA scale-free networks IpOl HH with the same size are 
first produced, and then these two networks are connected by 
few links. In order to normalize the terms of community net- 
work, we define the links between two communities as weak 
ties l27ll . and call the nodes connected by these weak ties brid- 
genesses 111 111 . Obviously, this network has a strong commu- 
nity structure because of few weak ties. With the increase of 
weak tie number, the community structure will be weakened. 

Network modularity Q, a popular evaluating indicator in 
measuring community structure 1 1], is defined as 



where l s and d s represent the number of intra-links and the 
sum of degrees of the nodes in community s respectively, L 
denotes the total link number in the network, and C is the 
number of communities. Here < Q < 1, the larger Q is, 
the stronger community structure is. However, this Q can not 
accurately characterize the community strength of the network 
with two communities ll3ltl . To address this shortcoming, the 
normalized Q n is defined as 
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where Q ra nd corresponds to random network with the same 
degree sequence, and Q ma x is the modularity of the network 
without inter-community links, i. e., Iab = 0. After this 
normalization, Q n is range from to 1. 



B. Dynamic Process 

Owing to the simplicity of SI model, the effects of dif- 
ferent contact patterns on epidemic spreading can be clearly 
understood. Although other disease models such as SIS are 
even more practical, more parameters such as the recovered 
rate /x in SIS model make epidemic dynamics more compli- 
cated. Take the threshold of outbreaks for example. Different 
initial seeds can result in distinct thresholds due to their lo- 
cal structures such as degree. For simplicity, we only study 
susceptible-infected (SI) 111411 spreading dynamics in contact 
process (CP) ll52l - l62ll through numerical simulations. In SI 
model, '5" and T represent respectively the susceptible (or 
healthy) state and the infected state. At the beginning, a node 
is selected as the initial infected (i. e., seed) and all other 
nodes are in '5' state. At each time step, each infected node 
randomly contacts one of its neighbors, and then the contacted 
neighboring node will be infected with probability A if it is in 
the healthy state, or else it will retain its state. Once an indi- 
vidual is infected, it will keep its state forever. To eliminate 
the stochastic effect of disease propagation, we set A = 1. 



C. Statistical Parameter 

In view of the relative independence of a local community, 
we take the dynamics and its variability into account. When 
a disease emerges in community network, it is very important 
for a local community to keep a watchful eye on two statisti- 
cal parameters: the arrival time and the prevalence of disease, 
where the arrival time of disease t a is defined as the moment 
that infectious individual first occurs in the community in each 
realization, and the prevalence i(t) is the density of infected 
individuals at time t. In order to investigate the predictability 
of epidemic dynamics, the variability of arrival time (preva- 
lence) is defined as the relative variation of the arrival time 
(prevalence) given by (HSU 5! 
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and 



A [i(t)] 



(<(*)> 



(4) 



A(t a ) = 0(A[z(t)] = 0) denotes all independent dy- 
namics realizations are essentially the same, and the arrival 
time (prevalence) in the network is deterministic. Larger 
A(f a )(A[i(f)]) means worse predictability that a particular 
realization is far from average over all independent realiza- 
tions. 



III. THE EFFECT OF WEAK TIES 



A. The Effect of Weak Tie with Different Degrees 

In consideration of the degree heterogeneity in this com- 
munity network, there may be different kinds of weak ties, 
that is, a pair of bridgenesses (i. e., 6^ and &#) connected 
by a weak tie may have different degrees k A and k B . In CP, 
weak ties with the different combinations of k A O k B have 
different effects on the propagation from the community A 
to the community B. On the other hand, in a network with- 
out community structure, different seeds have little impact on 
epidemic spreading lT4lTl . In other words, different k B almost 
don't affect epidemic spreading in the second community B. 
Therefore, we investigate how different k A influence the pre- 
dictability of epidemic dynamics in the second community B. 
As a first step towards this, we only consider the case of one 
weak tie between the community A and B, in which a node 
with the fixed degree in the first community A is connected 
with a randomly chosen node in the second community B. 

Fig. [T] shows the mean arrival time (t a ) and its variabil- 
ity A(t a ) in the second community B when different kinds 
of weak ties are created. When the bridgeness b A in the first 
community A is chosen as seed, i. e., d = 0, the mean arrival 
time (t a ) is linear with the degree of brigeness k A . As the 
bridgeness of the second community b B is one neighbor of the 
brigeness b A , it will be infected with probability 1 jk A at each 
time step, which is obviously a Poisson process. So the mean 
arrival time (t a ) is equal to k A and its relative variation is 
a(t a )/(t a ) ~ 1. When disease seed is a node with one step to 
the bridgeness b A (i. e., d = 1), the mean arrival time (t a ) for 
different k A will increase generally compared with the case 
ofd = 0. For small k A , (t a ) is significantly greater than that 
of d = 0, e. _ 

k A = 5. But for large k A , the relative change of (t a ) is very 



(t a ){d = 1) w 15 » (t a ){d = 0) « 5 for 
h - o. iuf large k A 
little, e. g., (t a )(d = 1) |a 248 > (t a )(d = 0) » 240 for 
k A = 242. The reason is that the infection of the bridgeness 
b B can be divided into two processes: the bridgeness b A is first 
infected in t x « 10, e. g., (t a )(d = l)-(t a )(d = 0) 10 for 
k A = 5 and (t a )(d = 1) - (t a ) (d = 0) « 8 for k A = 242, and 
then the bridgeness b B is infected by b A in t% € {t\,t\ + k A \. 
With the further increasing of distance of seed to the bridge- 
ness b A (such as d — 3,4), the mean arrival times (t a ) are 
nearly the same for the fixed k A . It can be understood that 
owing to the finite size effect of network with average shortest 
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FIG. 1: (color online). The mean arrival time (t a ) and its variabil- 
ity A(t a ) as a function of the degree of the bridgeness k A where 
the "squares", "circles", "triangleups", "triangledowns", and "dia- 
monds" denote the cases of seeds with d = 0, 1, 2, 3, and 4, respec- 
tively, (a) (t a ) versus k A , (b) A(t a ) versus k A . The parameters are 
chosen as N A = N B = 0.5 x 10 4 , (k) = 10, A = 1. We perform 
the experiments in 10 2 different networks, each of which are tested 
in 10 3 independent realizations. 



path length (L) ~ 3.7, the bridgenee b A is infected till overall 
outbreak emerges in the first community A. 

In Fig.Q](b), the variabilities of arrival time A{t a ) for dif- 
ferent k A axe. shown. When d = 0, A(i a ) for different k A are 
approximately equal to 1 because these infections are Poisson 
processes. Interestingly, compared with the case of d = 0, 
A(t a ) for d = 1 decrease generally, that is to say the further 
distance to the initial seed can lead to the better predictability 
of arrival time. As mentioned above, the infection of the brid- 
geness b B has two processes, and thus its variability A(t a ) 
can be written as 



A(f„) = A[h+t 2 ], 



(5) 



Where t\ , t% denote the time durations of the first process and 
the second process, respectively. Substituting it into Eq. (O, 
we obtain 



A(t„) 



(h+h) ' 



(6) 



where D{t\ + t 2 ) = ((ti + <2) 2 ) — (ti + ■ Considering 
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the independence of these two processes, Eq. © is reduced to 



A(i«) 



y/D{h)+D{t 2 ) 
(tl) + <*2> ' 



(7) 



where D(ti) and Dfe) are the time variance of the first pro- 
cess and the second process, respectively. 

In the first process, there are two basic spreading pathways 
through which the bridgeness b A may be infected: the brid- 
geness may be infected directly by the initial seed (i. e., 
the neighboring node j of the bridgeness b A ) with probabil- 
ity the other neighboring nodes are more likely to in- 
fect b A when the overall outbreak occurs in the first commu- 
nity. Although the first route is a Poisson process, the vari- 
ability A(£i) in t% will be less than 1 due to the determinacy 
of the second pathway. In the second process, the variability 
A (£2) — 1 because the infection in is a Poisson process. 
As D(t) = [A(t)(t)} 2 , we have 



A(t„) 



(tl) + <*2> 



(8) 



Obviously, A(t a ) must be less than 1 when d = 1. As 
(£2) = k A increases with k A , A(t a ) will also increase ac- 
cording to Eq. (0, which is verified by the results in Fig. Q] 
(b). Especially for the very large k A , the variability is very 
close to 1, This means that although the large degree of the 
bridgeness can delay the mean arrival time of disease, it cause 
the worst predictability of arrival time. When d > 2, A(ii) 
of the first process will be more determined. Thus, A(t a ) for 
small k A will become small, e. g., A(t a ) s» 0.31(rf = 4) < 
A(t a ) w 0.52(d = 1) for k A = 5. The above results demon- 
strate that when the degree of bridgeness is small, the further 
distance of the initial seed to the bridgeness can result in bet- 
ter predictability of arrival time due to the determinacy of the 
first process, while the variability A(t a ) — > 1 is almost not 
affected for the very large k A because of the Poisson property 
of the second process. 

Next, we focus on the statistical parameter A[i(t)]. 
Ref. 114111 showed that the variability of prevalence in a local 
community is very large at the beginning of outbreaks, which 
makes the prediction of prevalence hard. For this reason, we 
pay attention to the variability of prevalence in the early stage 
of outbreaks. In Fig. |2] (a), the mean prevalence (i(T)) at 
T = 20 decreases with k A and d, which is a complete rever- 
sal of (t a ) trend in Fig. [U(a). But Fig. |2](b) shows that its 
variability increases with k A and d, which is in accordance 
with the trend of (t a ). On the one hand, as the very large vari- 
ability of prevalence in the early stage is originated from the 
uncertain arrival time of disease 04111 . it is difficult for large 
k A (corresponding to large (t a )) to make sure the arrival of 
disease within T = 20. This will result in the large variability 
of prevalence, e. g., A[i(T)} 5.40(d = 0) for k A = 242. 
On the other hand, the further distance of the initial seed to 
the bridgeness b A also makes disease more difficult to arrive 
at the bridgeness b B , and can thus cause larger A[i(T)], e. 
g., A[i(T)] » 5.60(d = 4) > A[i(T)} w 0.60(d = 0) for 
k A = 5. These results imply that both the large degree of 




FIG. 2: (color online). At T = 20, the mean prevalence (i(T)) 
and its variability A[i(T)] as a function of the degree of bridgeness 
k A where the "squares", "circles", "triangleups", "triangledowns", 
and "diamonds" denote the cases of seeds with d = 0,1, 2, 3, and 
4, respectively, (a) (i(T)} versus k A , (b) A[i(T)] versus k A . The 
parameters are chosen as Na = Nb = 0.5 x 10 4 , (k) = 10, A = 1. 
We perform the experiments in 10 2 different networks, each of which 
are tested in 10 3 independent realizations. 



brigeness and the further distance of the initial seed can make 
the prediction of prevalence very hard. 



B. The Effect of Different Initial Seeds when d = 1 

From Sec. A, we know the degree heterogeneity of the 
bridgeness b A has a significant impact on the predictability 
of epidemic dynamics. In this case, we investigate the ef- 
fects of different degrees of the initial seed on the variabil- 
ity of epidemic dynamics in the second community B when 
the distance between the initial seed and the bridgeness b A 
is d = 1. As shown in Figs. [3] and [4] when k A = 5 is small, 
A(t a ) and A[i(T)] are obviously affected by the degree of the 
seed, while there is almost no effect on the variability when 
k A = 242 is large. 

When k A = 5, the large degree of the initial seed will result 
in large (t a ) (Fig.[3](a)) and small A(t a ) (Fig. 01(c)). Owing 
to the finite contact ability of the initial seed with large degree, 
it will cost more time to infect bridgeness b A in the first pro- 
cess, that is large (t a ). In this process, the bridgeness b may 
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FIG. 3: (color online). When the distance between initial seeds and 
birdgeness d = 1, the mean arrival time (t a ) and its variability A(t a ) 
as a function of the degree of the initial seed, where (t a ) versus k for 
kg = 5 (a) and kg = 242 (b), A(t a ) versus fc for kg = 5 (c) and 
fcjf = 242 (d). The results are averaged over 10 2 x 10 3 independent 
realizations in 10 2 networks. 
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FIG. 4: (color online). When d = 1, the mean prevalence (i(T)) and 
its variability A[i(T)] at T = 20 as a function of the de gree of the 
initial seed, where (i(T)) versus fc for kg = 5 (a) and fcf = 242 (b), 
A[i(T)] versus fc for fcf = 5 (c) and fcf = 242 (d). The results are 
averaged over 10 2 x 10 s independent realizations in 10 2 networks. 



be infected through two basic pathways. For the initial seed 
with small degree, the bridgeness b A is infected directly with 
higher probability. Thus, the first process introduces larger 
A(£i ) because of the randomness of Poisson process, and thus 
A(£ a ) in the whole process increases according to Eq. (|8). In 
addition, Fig.|4](c) shows that the variability of prevalence in- 
creases with the degree of the initial seed, which is consistent 
with the trend of the mean arrival time in Fig.[3](a). We should 
note that degree of the initial seed has an opposite effect on the 
variabilities of arrival time and prevalence, which may bring 
about a great trouble for pandemic prevention and control. 

When k A = 242, the whole infection process is domi- 
nated by the second Poisson process (i. e., very large (£2) 
in Eq. (0). Therefore, the different seeds with d = 1 can't af- 
fect the variability of epidemic dynamics visibly (see Figs. [3] 
(b),(d) andg](b),(d)). However, the very large A(t a ) « 0.90 
and A[i(T)] ~ 2.50 mean it is difficult to accurately forecast 
epidemic spreading when the bridgeness degree is large. In 
order to ensure the universality of the above results, other k A 
are also used to simulate this process. As expected, all simu- 
lations reveal the same conclusion. 



IV. THE EFFECT OF WEAK TIE NUMBER 

In real community networks with modularity Q 6 
[0.3, 0.7] 111], there are many weak ties between communi- 
ties. In this section, we would like to understand the effects of 
weak tie number on the predictability of epidemic dynamics. 
To gain a clear idea of the relation between modularity Q and 



weak tie number, Eq. (Q~|) is expanded as 

g -(£)«. TO 

S=l 

where Iab represents the amount of weak ties between the 
community A and the community B. When the network is 
connected randomly, Iab — (dAd>B)/2L, thus Q ran d — > 0; 
When Iab = and dA = ds, Q reaches the maximum value, 
i. e., Qmax — 0.5. Therefore Q can only range from to 
0.5. Substituting Qrand — and Q max — 0.5 into Eq. (O, 
we have 

Q„ = 2Q. (10) 

After this standardization, Q n can range from to 1. By 
adding weak tie number Iab between two communities ran- 
domly, we can obtain the community networks with different 

Qn- 

Fig. [5] shows the case of three kinds of initial seed: brige- 
ness, random node, and hub. In Fig. 0(a), with the increase of 
Q n , the mean arrival time (t a ) will increase because fewer 
weak ties lengthen the distance between two communities. 
Especially when Q n > 0.9, the mean arrival time (t a ) will 
increase rapidly. Compared with the other two cases, the case 
of brigeness chosen as the initial seed has the shortest (t a ). 
The case of random node includes the cases of brigeness and 
non-bridgeness, so (t a ) for random seed must be longer than 
that for bridgeness seed. For the case of hub, it has the longest 
mean arrival time. As the nodes connected by weak ties are 
chosen randomly, the nodes with small degree will be more 
probably chosen as bridgenesses due to the degree hetero- 
geneity, while it is very difficult for hubs to be bridgenesses. 
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FIG. 5: (color online). The mean arrival time (t a ) and its variabil- 
ity A(t a ) as a function of the modularity Q n where the "squares", 
"circles", and "triangles" denote the cases of bridgeness seed, ran- 
dom seed, and hub seed, respectively, (a) (t a ) versus Q n , (b) 
A(t„) versus Q n . The parameters are chosen as Na = Nb = 
0.5 x 10 4 , (k) = 10, A = 1. We perform the experiments in 10 2 
different networks, each of which are tested in 10 3 independent real- 
izations. 
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FIG. 6: (color online). At T = 2, the mean prevalence (i(T)) and 
its variability A[i(T)] as a function of as a function of the modu- 
larity Q„ where the "squares", "circles", and "triangles" denote the 
cases of bridgeness seed, random seed, and hub seed, respectively, 
(a) (i(T)) versus Q n , (b) A[i(T)] versus Q n . The parameters are 
chosen as N A = N B = 0.5 x 10 4 , (k) = 10, A = 1. We perform 
the experiments in 10 2 different networks, each of which are tested 
in 10 3 independent realizations. 



When the degree of bridgeness is small, the initial seed with 
large degree leads to the longer (t a ) (see Fig.|3](a)). 

Moreover, the variabilities of arrival time A(t a ) for these 
three cases are compared in Fig.|5](b). As mentioned in Fig.Q] 
(b), the further distance of the initial seed to the bridgeness 
b A can reduce the variability A(t a ). With the increase of Q n , 
the distance between two communities is lengthened by fewer 
weak ties, and thus A(t a ) will decrease. For example, when 
Q n increase from 0.9 to 1, A(t a ) for the random case decrease 
from 0.33 to 0.18 rapidly. For the case of hub, its arrival time 
has the most accurate predictability. This is because that the 
initial seed with large degree can lead to the low variability 
when the brigeness b A has small degree (see Fig.|3](c)). More 
significant, A(t a ) for the case of brigeness increases with the 
network modularity Q n , which is opposite to the other two 
cases. Even though a bridgeness bf is chosen as the initial 
seed, the infection of community B is not always through the 
weak tie of bridgeness bf because of the existence of many 
weak ties. In this case, there are two optional spreading path- 
ways towards the community B: the weak tie of bf and the 
other weak ties. From this viewpoint, the actual path length 
of epidemic spreading must be greater than 1. In other words, 
more links will result in the longer distance of the infection 
process. Therefore, more weak ties can reduce A(t a ) due to 
the increase of distance between the community A and the 
community B. Actually, more weak ties can increase the de- 
terministic of the second optional pathway, and thus enhance 
the predictability of arrival time. 

Furthermore, the effects of weak tie number on the pre- 
dictability of prevalence in the early stage of outbreaks are 
also analyzed in Fig. [6] As the mean arrival time increases 
with Q n in Fig. [5] (a), the prevalence (i(T)) at T = 2 
will decrease accordingly in Fig. [6] (a). (i(T)){hub) < 
(i(T)) (random) < (i(T)) (bridgeness) is resulted from 



(t a ) (hub) > (t a )(random) > (t a ) (bridgeness). Fig.|6](b) 
shows that the variability of prevalence A[i(T)} increases with 
Q n . For the case of bridgeness, the change of A[i(T)] is very 
little, which means the bridgeness palys a significant role in 
enhancing the predictability of prevalence 114 ill . For the cases 
of randomly chosen node and hub, A[i(T)] increases slowly 
when Q n E [0,0.9), while A[i(T)] increase rapidly when 
Q n G [0.9, 1) (e. g., A[i(T)] » 13.10 for Q n « 0.99), which 
is in accordance with the trend of (t a ). The results at other T 
value (e. g., T = 5, 10) reveal the same conclusion. It implies 
that strong community structure can increase the difficulty of 
the predictability of prevalence. 



V. CONCLUSIONS 

In conclusions, we have studied the effects of weak ties on 
the predictability of epidemic dynamics in the local commu- 
nity with degree heterogeneity. First, we have shown that the 
degree of bridgeness can remarkably influence the variabil- 
ities of both the arrival time and the prevalence of disease. 
With the increase of the degree of bridgeness, the mean ar- 
rival time and the outbreak (i. e., prevalence) of disease will 
be delayed, but their variabilities will also increase. In addi- 
tion, we have also shown that the distance of the initial seed 
to the bridgeness has different impacts on the epidemic pre- 
dictability under different conditions. When the degree of the 
bridgeness is small, the further distance of the initial seed to 
the bridgeness will enhance the predictability of arrival time, 
while the predictability of prevalence in the early stage will 
get worse. When the degree of the bridgeness is large, the 
variability of arrival time is almost close to 1 because of the 
Poisson property of the infection process, while the variability 
of prevalence is very large due to the uncertain arrival time of 
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disease. Second, we have investigated the effects of the initial 
seeds with different degrees on the variability of epidemic dy- 
namics when the distance between the initial seed and bridge- 
ness is equal to 1. When the degree of the bridgeness is small, 
the large degree of the initial seed will enhance the predictabil- 
ity of arrival time, while the predictability of prevalence in the 
early stage will be worse. When the degree of the bridgeness 
is large, the variability of epidemic sprading is almost not af- 
fected. Moreover, we also have analyzed the effects of weak 
tie number on the epidemic predictability where the results of 
three different initial seeds (i. e., brigeness, random node, and 
hub) were compared. With the increase of the network modu- 
larity, which is caused by the decrease of weak tie number, the 
variabilities of arrival time for the case of random node and 
hub will first decrease slowly, and then decrease rapidly as 
the community strength is very strong. It's important to note 
that the variability of prevalence will increase rapidly when 
the community strength becomes stronger, which is contrast 
to the trend of the variability of arrival time. By contrast, the 



variability of prevalence for the case of bridgeness increases 
with the network modularity, which displays the worst pre- 
dictability of arrival time. However, the best predictability of 
prevalence is observed when the bridgeness is first infected. 
The above results show that the variability of arrival time al- 
ways marks a complete reversal trend of that of prevalence, 
which implies it is impossible to predict epidemic spreading 
in the early stage of outbreaks accurately. 



Acknowledgments 

This work is supported by the NNSF of China (Grants Nos. 
11105025, 90924011), China Postdoctoral Science Founda- 
tion (Grant No. 20110491705), the Specialized Research 
Fund for the Doctoral Program of Higher Education (Grant 
No. 201 10185120021), and the Fundamental Research Funds 
for the Central Universities (Grant No. ZYGX201 1 J056). 



[1] M. E. J. Newman, Phys. Rev. Lett. 89, 208701 (2002); Phys. 

Rev. E 67, 026126 (2003); M. E. J. Newman and J. Park, Phys. 

Rev. E 68, 036122 (2003); M. E. J. Newman and M. Girvan, 

Phys. Rev. E 69, 0261 13 (2004). 
[2] Facebook, http://www.facebook.com 
[3] Youtube, http://www.youtube.com 
[4] Xiaonei, http://www.xiaonei.com 

[5] H. -H. Jo, R. K. Pan, and K. Kaski, PLoS ONE 6(8), e22687 
(2011). 

[6] E. Ferrara, P. De Meo, G. Fiumara, and A. Provetti, Procedia 

Comput. Sci. 00, 1 (2012). 
[7] P. A. Grabowicz, J. J. Ramasco, E. Moro, J. M. Pujol, and V. 

M. Eguiluz, PLoS ONE 7(1), e29358 (20 12). 
[8] M. S. Granovetter, Am. J. Sociol. 78, 1360 (1973). 
[9] J. Scott, Sociology 22, 109 (1988). 
[10] E. L. Berlow, Nature 398, 330 (1999). 

[11] X. -Q. Cheng, F. -X. Ren, H. -W. Shen, Z. -K. Zhang, and T. 
Zhou, J. Stat. Mech. P10011 (2010). 

[12] P. Csermely, Weak links: The universal key to the stability of 
networks and complex systems, (Heidelberg: Springer, Ger- 
many, 2009). 

[13] N. T. J. Bailey, The Mathematical Theory of Infectious Dis- 
eases, second ed., (Griffin, London, 1975). 

[14] R. M. Anderson and R. M. May, Infectious Disease of Humans, 
(Oxford University Press, Oxford, 1992). 

[15] O. Diekmann and J. A. P. Heesterbeek, Mathematical Epidemi- 
ology of Infectious Diseases: Model Building, Analysis and In- 
terpretation, (Wiley, New York, 2000). 

[16] D. J. Dailey and J. Gani, Epidemic Modelling: An Introduction, 
(Cambridge University Press, Cambridge, 2001). 

[17] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, D.- U. Hwang, 
Phys. Rep. 424, 175 (2006). 

[18] S.N. Dorogovtsev, A. V. Goltsev AV, J. F. F. Mendes, Rev. Mod. 
Phys. 80, 1275 (2008). 

[19] A. Barrat, M. Barthelmy, and A. Vespignani, Dynamical Pro- 
cesses on Complex Networks (Cambridge University Press, 
New York, 2008). 

[20] M. E. J. Newman, Networks: An Introduction (Oxford Univer- 
sity Press, Oxford, U.K., 2010). 



[21] M. Barthelemy, Phys. Rep. 499, 1 (2011). 
[22] Z. Liu and B. Hu, Europhys. Lett. 72(2), 315 (2005). 
[23] L. Huang, K. Park, and Y. -C. Lai, Phys. Rev. E 73, 035103(R) 
(2006). 

[24] Y. Zhou, Z. Liu, and J. Zhou, Chin. Phys. Lett. 24, 581 (2007). 
[25] G. Yan, Z.-Q. Fu, J. Ren, and W.-X. Wang, Phys. Rev. E 75, 

016108 (2007). 
[26] X. Wu and Z. Liu, Physica A 387, 623 (2008). 
[27] J. -P. Onnela, J. Saramaki, J. Hyvonen, G. Szabo, D. Lazer, K. 

Kaski, J. Kertesz, and A. -L. Barabasi, Proc. Nat. Acad. Sci. 

104, 7332 (2007). 
[28] D. Centola and M. Macy, Am. J. Sociol. 113, 702 (2007). 
[29] J. C. Zhao, J. J. Wu, andK. Xu, Phys. Rev. E 82, 016105 (2010). 
[30] A. Mihalik and P. Csermely, PLoS Comput. Biol. 7(10), 

el002187 (2011). 

[31] H. Yang, M. Tang, H. -F. Zhang, |arXiy:1205.4352V l (2012). 
[32] L. Hufnagel, D. Brockmann, and T. Geisel, Proc. Natl. Acad. 

Sci. U. S. A. 101, 15124 (2004). 
[33] M. Barthelemy, A. Barrat, R. Pastor-Satorras, and A. Vespig- 
nani, J. Theor. Biol. 235, 275 (2005). 
[34] A. Gautreau, A. Barrat, and M. Barthelemy, J. Stat. Mech. 2007, 

L09001. 

[35] A. Gautreau, A. Barrat, and M. Barthelemy, J. Theor. Biol. 251, 
509 (2008). 

[36] M. Barthelemy, C. Godreche, and J. -M. Luck, J. Theor. Biol. 
267, 554 (2010). 

[37] V. Colizza, A. Barrat, M. Barthlemy, and A. Vespignani, Proc. 

Natl. Acad. Sci. U. S. A. 103, 2015 (2006). 
[38] P. Crepey, F. P. Alvarez, and M. Barthelemy, Phys. Rev. E 73, 

046131 (2006). 

[39] M. Loecher and J. Kadtke, Phys. Lett. A 366, 535 (2007). 

[40] W. Huang and C. Li, J. Stat. Mech. 2007, P01014. 

[41] K. Gong, M. Tang, H. Yang, and M. Shang, Chaos 21, 043130 

(2011). 

[42] Z. -D. Zhao, Y. Liu, and M. Tang, arXiv: 1204.0706V 1 (2012). 
[43] J. D. Han, N. Bertin, T. Hao, D. S. Goldberg, G. F. Berriz, L. V. 

Zhang, D. Dupuy, A. J. Walhout, M. E. Cusick, F. P. Roth, and 

M. Vidal, Nature 430, 88 (2004). 
[44] N. N. Batada, T. Reguly, A. Breitkreutz, L. Boucher, B. -J. Bre- 



8 



itkreutz, L. D. Hurst, and M. Tyers, PLoS Biol. 4, e317 (2006). 
[45] N. N. Batada, T. Reguly, A. Breitkreutz, L. Boucher, B. -J. Bre- 

itkreutz, L. D. Hurst, and M. Tyers, PLoS Biol. 5, el54 (2007). 
[46] Y. Wang, Y. -Q. Hu, Z. -R. Di, and Y. Fan, Physica A 390, 4027 

(2011). 

[47] M. Girvan and M. E. J. Newman, Proc. Natl. Acad. Sci. U. S. 

A. 99, 7821 (2002). 
[48] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D.-U. 

Hwang, Phys. Rep. 424, 175 (2006). 
[49] M. E. J. Newman, Nat. Phys. 8, 25 (2012). 
[50] A. -L. Barabasi and R. Albert, Science 286, 509 (1999). 
[51] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 (2000). 
[52] T. Zhou, J. -G. Liu, W. -J. Bai, G. -R. Chen, and B.-H. Wang, 

Phys. Rev. E 74, 056109 (2006). 
[53] C. Castellano and R. Pastor-Satorras, Phys. Rev. Lett. 96, 

038701 (2006). 

[54] M. Ha, H. Hong, and H. Park, Phys. Rev. Lett. 98, 029801 



(2007). 

[55] C. Castellano and R. Pastor-Satorras, Phys. Rev. Lett. 98, 
029802 (2007). 

[56] H. Hong, M. Ha, and H. Park, Phys. Rev. Lett. 98, 258701 

(2007) . 

[57] R. Yang, B. -H. Wang, J. Ren, W. -J. Bai, Z. -W. Shi, W. -X. 

Wang, T. Zhou, Phys. Lett. A 364, 189 (2007). 
[58] R. Yang, T. Zhou, Y. -B. Xie, Y. -C. Lai, and B. -H. Wang, Phys. 

Rev. E 78, 066109 (2008). 
[59] R. Yang, L. Huang, and Y. -C. Lai, Phys. Rev. E 78, 026111 

(2008) . 

[60] J. D. Noh and H. Park, Phys. Rev. E 79, 0561 15 (2009). 

[61] S. H. Lee, M. Ha, H. Jeong, J. D. Noh, and H. Park, Phys. Rev. 

E80, 051127 (2009). 
[62] M. A. Munoz, R. Juhasz, C. Castellano, and G. Odor, Phys. 

Rev. Lett. 105, 128701 (2010). 



